ZEROS OF PARTIAL SUMS OF POWER SERIES FOR A FAMILY OF 
EXPONENTIAL INTEGRALS, I 



ANTONIO R. VARGAS 



Abstract. We are interested in studying the asymptotic behavior of the zeros of partial 
sums of power series for a family of entire functions defined by exponential integrals. The 
zeros grow on the order of O(n), and after rescaling we explicitly calculate their limit 
curve. We find that the rate that the zeros approach the curve depends on the order of 
the singularities/zeros of the integrand in the exponential integrals. As an application 
of our findings we derive results concerning the zeros of partial sums of power series for 
Bessel functions of the first kind. 



1. Introduction 

In this paper we are concerned with the asymptotic behavior of the zeros of the polynomial 
sequence given by the partial sums of a convergent power series. If / is a function which is 
analytic at the origin, then it can be represented by a power series 

oo 
k=0 

which converges near z — 0. We denote the n th partial sum of this power series by 

n 

(1) s n [f](z) = ^2a k z k , 

k=0 

and we refer to s n [f](z) as the n th section of /. 

If the power series for / has a finite radius of convergence, then it is a classical result 
of Jentzsch [8] that every point on the circle of convergence of the power series will be a 
limit point of the zeros of the sections s n [f](z). In addition, because power series converge 
uniformly on compact subsets of their domains of convergence, Hurwitz's theorem (see, e.g., 
[9, p. 4]) tells us that any zero of / inside the radius of convergence will also be a limit point 
of the zeros of the sections. 

The behavior of the zeros becomes much more interesting when / is entire. The topic 
can be traced back to Szego, who studied in [14] the sections of the exponential function e 2 , 
given by 

n 

s n [exp](» = tt • 

k=0 

Szego found that the zeros of the normalized sections s„[exp](nz) have as their set of limit 
points the simple closed loop 

D = {z G C : \z\ < 1 and l^e 1 ^ = l} . 

This is often referred to as the Szego curve. The rate that the zeros approach this curve 
was first studied by Buckholtz [5], who showed that every zero of s„ [exp] (nz) lies within a 
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distance of 2e/^/n of D. Carpenter, Varga, and Waldvogel [6] examined this phenomenon 
in detail and showed that Buckholtz's result gives the best-possible asymptotic order. The 
statement of the theorem involves the complementary error function 

2 f°° _ f 2 

erfc z = —= I e dt, 

where the path of integration begins at z and travels to the right to oo. 

Theorem 1.1 (CVW). For fi C C, define maxdist(ft, C) = sup 2GO {dist(z, C)}. // 
{ z k,n}k=i are the zeros of s„[exp](nz) and if t\ is the zero of the complementary error 
function erfc closest to the origin in the upper half-plane, then 

liminf v 7 ^- maxdist({z fc .„}^ =1 ,D) > Ufa) + $sfa) ~ 0.636657. 

n— >oo 

The authors also showed that the zeros which are bounded away from the point z = 1 
approach the Szego curve more quickly. 

Theorem 1.2 (CVW). Let Cs be the ball of radius S centered at z = 1. If {zk, n }k=i are 
the zeros of s n [exp] (nz) and if S is any fixed number with < S < 1, then 

maxdist({z fei „}fc =1 \ C s , D) = O 

as n — > oo. 

There have been some similar results for other power series, notably sine and cosine (see 
[15] and the references therein), a class of confluent hypergeometric functions [11], finite 
sums of exponentials [2] , and even some classes of divergent power series [7] . For a survey 
of this topic the reader is referred to [16]. 

In this paper we will continue in this vein, studying power series defined by exponential 
integrals of a certain form, which we introduce in Section 2. Special cases of these functions 
include the aforementioned confluent hypergeometric functions as well as Bessel functions 
of the first kind, the latter of which we will discuss in detail in Section 4. 

Our results were particularly inspired by the work of Norfolk [11] on the confluent hy- 
pergeometric functions, defined by 

^ftM = r(6)gj^. 

Norfolk studied the case where b is real and 6^1,0, —1, —2, . . . (the function is defined by 
analytic continuation when b < 0). His main tool was an exponential integral representation 
of the functions valid for b > 1 (see Section 4). The extension of some of Norfolk's results 
to the case where b is complex with U(b) > 1 will be a consequence of our main theorems. 

This paper is the first in a series of two. The results within were originally obtained in 
the author's Master's thesis [16] and form an analogue to Theorem 1.2. The second paper 
in the series will contain new results; in particular an analogue to Theorem 1.1. 

2. Definitions and preliminaries 

The functions we are interested in are defined by integrals of the form 
ip(t)e zt dt. The restrictions we place on the function ip are determined essentially by 
the abilities of Watson's lemma, discussed below. 
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Suppose < a, b < oo and let ip : [—a, b] — > C U {00} be a function satisfying 

rb 

\tp(t)\dt < 00 

and ip(t) = {t + aYfi{t + a) = (b - t) v f 2 (b - t), where 

(1) (i,i/6C with > -1 and 5ft(i/) > -1, 

(2) /1, /2 : [0, a + b] — >• C U {00} with /i(0) and /2(C)) both finite and nonzero, 

(3) in a neighborhood of t = 0, both /((i) and /^(t) exist and are bounded. 
Define 

F(jb) = f p(t)e zt dt. 

J —a 

It is a consequence of the dominated convergence theorem that the function F is entire, and 
its sections are given by the formula 

™ k rb 

(2) »»[*!(*) = £w / 

We can view -F(z) as the exponential generating function of these particular integral mo- 
ments of If. 

Properties (1), (2), and (3) above describe how the function tp behaves near the endpoints 
of integration, and perhaps more importantly they name various quantities we will refer to 
throughout the paper. Property (3) serves a special purpose: it allows us to determine a 
simple error term in the asymptotic expansion of F through the use of Watson's lemma 
(Theorem 2.2). 

Next we collect some theorems which will aid us in proving our results. The first such 
theorem is due to Rosenbloom [12, 13]. 

A sequence of sections {sat[/](z)} is said to have a positive fraction of zeros in any sector 
with vertex at the origin if 

liminf M%M >0 

for any fixed 6\ and 62, where ^(^i, #2) is the number of zeros of the section s?f[f](z) in 
the sector Q\ < argz < 62 ■ Roscnbloom's result is as follows. 



Theorem 2.1 (Rosenbloom). Let 



k=0 



be an entire function of finite positive order p and let s n [f](z) be as in equation (1). Define 
jj°(i?) to be the number of zeros of s n [f](z) in the disk \z\ < R and define p n = |a n | -1 /". 
For any increasing sequence of indices {N} such that the sequence of sections {sN[f](z)} 
has a positive fraction of zeros in any sector with vertex at the origin, the existence of which 
is guaranteed, and for < r < 1 and e > we have 

liminf ^ '- > 1 - r p > 0. 

N^oo N 

Further, the number of zeros o/sat[/](z) satisfying 

\p N z\ > e 1/p + e 

is bounded. 
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We will also require Watson's lemma on the asymptotic behavior of exponential integrals. 
We refer the reader to [10] for a thorough discussion of this result. In the following, A is a 
complex parameter. 

Theorem 2.2 (Watson's lemma). Suppose < T < oo and (p : [0,T] — > C U {oo} is a 
function satisfying 



Jo 

and ip(t) — t a h(t), where 5ft(er) > —1, h(0) ^ 0, and h'(t) exists and is bounded in a 
neighborhood oft = 0. Then the exponential integral 



as A — > oo with | arg A| < 9 for any fixed < 9 < ir/2. 

Though this form of Watson's lemma only gives an asymptotic for $(A) as A — > oo to the 
right, it can easily be extended to address the case when A — > oo to the left if we assume 
that T is finite. 

Corollary 2.3. Suppose < T < oo and ip : [0,T] — > C U {oo} is a function satisfying 



Jo 

and ip(t) = (T — t) a h(T — t), where 5R(ct) > —1, h(0) ^ 0, and h'(t) exists and is bounded 
in a neighborhood oft = 0. Then the exponential integral 





is finite for all 5R(A) > ; and 



$(A) 



fo(0)I> + l) 
A CT +! 



+ 0(A- ff - 2 ) 





is finite for all 5R(A) > ; and 




as A — > oo with \ arg A| < 8 for any fixed < 9 < ir/2. 



Proof. We have 




as A — > oo with | arg A| < 9 for any fixed < 9 < n/2, by Watson's lemma. 



□ 
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3. Main results 
Recall from Section 2 that we are concerned with functions of the form 

rb 



F(z) = f <p(]t), 

J —a 



zt 



dt 



with ip satisfying some light requirements. The n section of F is the polynomial 

rb 

tp{t)t k dt. 



n k rb 
fc=0 J ~ a 



The statements and proofs of the main results of this section depend on the relative sizes 
of a and b and of 3?(/x) and To this end, define 

c = max{a, b} 

and 

'JR(/Lt) if a > b, 

SR(i/) if a < b, 

min{sft(//),sft(i/)} if a = b. 

Lastly, let {N} be an increasing subsequence of the indices {n} for which the sequence of 
sections {sn[F](z)} has a positive fraction of zeros in any sector with vertex at the origin 
(see Section 2). If a = b and 9?(/i) = 3f?(^) , we also impose the condition that the indices 
{N} are chosen so that quantity 

(3) (-1)^/1 (0)r(/i + 1) + / 2 (0)I> + l)a»-»N^ 

is bounded away from 0. Such a subsequence is guaranteed to exist by Theorem 2.1. This 
condition imposed on the quantity in (3) ensures that we can use the asymptotic represen- 
tations derived in Lemma 3.5 and Lemma 3.6 without incident. 

Theorem 3.1. Let F be an exponential integral function as in Section 2 and let {N} be a 
subsequence of the indices {n} as defined above. 

(i) Let {zjv} be a sequence of complex numbers such that sn[F](Nzn) = for all N 
which has a limit point in the region —1/a < Re(z) < 0. Then the elements of the 
sequence satisfy 



1+aZN I = i + (e - »o*) + l) l ^r + 0{1/N) 



\czNe 
as N — > oo. 

(ii) Let {zn} be a sequence of complex numbers such that sn[F](Nzn) = for all N 
which has a limit point in the region < 5R(z) < 1/6. Then the elements of the 
sequence satisfy 

Iczve 1 -^ | = 1 + (e - »M + l)^JT + °( W 

as N — > oo. 
Theorem 3.2. It is true that 
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Z < 



1 



(i) Every point on the curve 

D a , b = \z e C : $t(z) < 

J < z £ C : > 0, \z\ < 



and \cze 
1 

ana 



l+az\ 



= 1 



is a limit point of the zeros of the sections s n [F] (Nz) . 
(ii) The limit points of the zeros of the sections sj\r [F] {Nz) on the imaginary axis are 
precisely the points on the line segment 



D 



imag 



= <zeC; 5R(z) = and \z\ 



< 



1 



(Hi) All but finitely many limit points of the zeros of the sections sn[F](Nz) lie on 
D at bUDi ma g. Any limit points not on D a ^UDi mag lie outside of the circle \z\ = 1/c. 

Remark 3.3. The expressions in parts (i) and (ii) of Theorem 3.1 give information about 
whether the zeros eventually lie on the inside or the outside of the limit curve based on signs 
of the quantities £-&(/z) + l/2 and £ - ft(iy) + 1/2. For example, if + 1/2 > 0, then 

it will eventually be true that the zeros which approach D a (, in the left half-plane will satisfy 
\cze 1+az | > 1 and hence will lie outside of -D ,b- If either of the quantities £ — + 1/2 or 
£ — 5i(V) + 1/2 is zero then the theorem does not give any information about the direction 
from which the zeros approach the relevant part of the curve. 

Figure 1 and Figure 2 illustrate two interesting cases of the results in Theorem 3.1 and 
Theorem 3.2. The reader may refer to [16] for more plots of this kind. 



0.4 



0.2 



0.0 



-0.2 



-0.4 




-1.0 -0.5 0.0 0.5 1.0 



Figure 1. Zeros of the normalized section sso [F](80z) and the limit curve 
D a ,b U Amag with a = b = 1 and ip(t) = (1 - tf/ 2 (l + t)- 1 / 2 ^. Note 
that £ — $l(v) + 1/2 = —3/2 < 0, which predicts that the zeros in the right 
half-plane will approach the limit curve from the interior. 



To prove these theorems we will require a few lemmas, the first of which concerns the 
asymptotic behavior of the function F. 
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-0.5 0.0 0.5 1.0 



Figure 2. Zeros of the normalized section sgo[i ;l ](80z) and the limit curve 
Ax.fcUAmag with a = 1, b = §, and tp(t) = (§ - + i) 1/2 ~ 4 . Note that 
£ — $t(fi) + 1/2 = 0, which predicts that the zeros in the left half-plane will 
approach the limit curve at a rate of 0(l/n). 

Lemma 3.4. As n — > oo, 

F(nz) = / 1 (0)r(/i + l)(-ri2)-' t - 1 e- onz (l + O(l/n)) 
when z is restricted to a compact subset of 3?(z) < 7 and 

F(nz) = f 2 (0)T(u + l)(nz)- u - 1 e bnz (l + 0(l/n)) 
when z is restricted to a compact subset o/5R(z) > 0. 

Proof. This follows from a direct application of Corollary 2.3. To see this, suppose first that 
z is restricted to a compact subset of 3t(z) < 0, and make the substitution t = b — s in the 
integral for F(nz) to get 

/b f-a+b 
Lp{t)e nzt dt = e bnz / ip{b - s)e- nzs ds, 
-a JO 

which, after replacing z with — z, is of the form required by the corollary. Next, suppose 
that z is restricted to a compact subset of 5ft(z) > 0, and make the substitution t = s — a in 
the definition of F(nz) to get 

/b ra-\-b 
tp{t)e nzt dt = e~ anz / <p(s - a)e nzs ds, 
-a JO 

which is also of the required form. □ 

We must also find asymptotics for the integral moments of tp and hence for the power 
series coefficients of F. 
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Lemma 3.5. We have 

<p(t)t n dt = (-l) n /i (O)r(Ai + l)n-"- 1 o n+ ' ,+1 + 0(n->*- 2 a n ) 



+ / 2 (0)r(i/ + l)n- / - 1 b n+l,+1 + 0(n- v - 2 b n ) . 



as n — > oo. 



Proof. If a 7^ we calculate 

tp{t)t n dt = {-a) n I ip{t)e nloe( ~ t/a) dt 



(-a) n / <p(s-a)e nlog(1 - s/a) ds 
Jo 

(-a) n [ s' i / 1 (s)e" lo6 ( 1 - s / a ^s. 



Letting s = a(l — e ' ) gives 

<p(t)t n dt = (-l)" a n+ " +1 / (1 - e-'y/i(a - ae - r )e- r e- nr dr 



= (-l) n a n+ >* +1 / r"ip a (r)t 
Jo 



where 

Mr) = " ae- r )e- r 

has a bounded derivative in a neighborhood of r = 0. We may now apply Watson's lemma 
to conclude that 

ip{t)t n dt = (-l) n %(0)T(fi + l)n-^- 1 a n+ti+1 + 0(n-^- 2 a n ) 

= (-l)"/l(0)r(/i + l)n-^- l a n+ ^ +1 + 0(n-^-' 2 a n ) . 

Using an identical argument we find that 
i 

(p(t)t n dt = / a (0)r(i/ + l)n-"-' L b n+v+l + 0(n- u - 2 b n ) , 

which completes the proof. □ 
A similar argument can be used to prove the following. 



Lemma 3.6. 

rb 



A- t n+1 dt = /l( ° )r(Af + 1) n-"- 1 a n+ ^+ 2 + 0(n-»- 2 a n ) 

1 - zt K ' 1 + az v ' 



n 



l-bz 

as n — > oo uniformly when z is restricted to a compact subset of the doubly-punctured plane 
{z e C : z ^ -1/a and z ^ l/b}. 

In the next lemma we will use the result of Lemma 3.5 to place a bound on where the 
limit points of the zeros of the normalized sections sn[F](Nz) may lie. 
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Lemma 3.7. All but finitely many limit points of the zeros of the sections s^[F](Nz) lie 
in the disk \z\ < 1/c, where c — max{a, b}. 

Proof. From Stirling's formula 







we have 



. e 
(n!)V 



as n ~ > oo, and with the aid of Lemma 3.5 we calculate 



<p(t)r dt 



-l/N 



as N —> oo, where the subsequence of indices {N} is as defined above Theorem 3.1. Com- 
bining these we see that the power series coefficients of F(z), which are given by 

rb 



satisfy 



1 



Pn = |ojv| 



(p(t)t n dt, 



-l/N 



N 



where pjv is as defined in Theorem 2.1. Thus the order p of F is calculated to be 

p = lim sup = 1 , 

n^oo fog p n 

where all indices n are taken into account (see, e.g., [3, p. 9]). 
By Theorem 2.1 we have, for any e > 0, 



lim inf 

7V->oo 



tt^((l + e)JV/c) _ 
N 



where ttjv((l + e )N/c) is the number of zeros of Sn[F](z) in the disk \z\ < (1 + e)N/c. 
From this we conclude that all but finitely many limit points of the zeros of the sections 
sn[F](Nz) lie in the disk \z\ < 1/c, as desired. □ 

We may now prove the first theorem in this section. 

Proof of Theorem 3.1. By definition we have 



and 



F(nz) 
s n [F](nz) = 



Subtracting these we get 

F(nz)-s n [F](nz) 



<p(t)e nzt dt, 



(p(t)s n [exp](nzt) dt. 



ip(t) (e nzt - s n [exp](nzt)) dt 



(4) 

where 



— a 
b 



<p(t)e nzt g n (zt) dt, 
g n (z) = 1 - e~ nz s n [ex.p](nz). 
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It was shown by Szego in [14] (see also [6], [4], and [11]) that 

i l— z\ n 

(5) 9n(z) = ^=^- ■ ^(l- tn(z)), 

V27rn 1 - z \ / 

where e n (z) — 0(l/n) as n — > oo uniformly when z is restricted to a compact subset of 
3?(z) < 1. Upon substituting this into equation (4) we get 

F(nz) - s n [F}(nz) = ^= [ A- t n+1 (l - e n (zt)) dt. 

\/2nn J - a 1 — Zt \ / 

If c = max{et, b} it follows that zeros of s n [F\(nz) which remain in compact subsets of the 
region \3t(z)\ < l/c satisfy 

(6) F(nz) = f A f + 1 dt (l + 0(l/n)) 

asn-> oo, where we have used Lemma 3.6 to bring the error term outside of the integral. 

Let {N} be a sequence of indices as defined above Theorem 3.1. We showed in Lemma 3.7 
that all but finitely many the limit points of the zeros of the sections s n [F] (Nz) lie in the 
disk \z\ < l/c. Define 

A = \z € C : \z\ < -, z^--, and 

[ cab 

Suppose first that {z^} is a sequence in C such that sn[F](Nzn) = for all TV and such 
that the sequence has a limit point in A n {z £ C : 5R(z) < 0}. Note that this implies there 
is a 6 > such that \zn + l/ a l > S for N large enough. 

It follows from Lemma 3.4 that 

(7) \F(Nz N )\V N = |e—*| (l - (»(„) + 1) + 0(l/iV)) 
as TV — > oo, and if 

{5K(/i) if a > b, 

if a < 6, 
nun{3ft(/i),5ft(i/)} if a = 6, 



then wc have from Lemma 3.6 that 



(8) 



. t N+1 dt 

\/2ttN J - a 1- zt 



= <<--,\il-i^+l) l °jf+0(l/N) 



as TV — > oo. Upon substituting equations (7) and (8) into equation (6) we see that these 
zeros zjv satisfy 



\CZ]ye 



1+aZN I = i + (e - »o*) + + 0(1/A ° 



as TV — > 00, which proves part (i) of Theorem 3.1. 

Suppose now that {zjy} is a sequence such that Sn[F](Nzn) — for all N and such that 
the sequence has a limit point in A n {z € C : > 0}. Note that this implies there is a 
S > such that |zjv — 1/6| > # for N large enough. 

Here it follows from Lemma 3.4 that 

\F(Nz N )\^ N = \e b ^\ (l - 0R(t/) + 1) ^ + 0(1/JV) 
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as N — > oo. Subtituting this and equation (8) into equation (6) we see that these zeros zn 
satisfy 

Ic^e 1 '^ | = 1 + U ~ «M + 1)^ + °^I N ) 

as N — > oo, which proves part (ii) of Theorem 3.1. □ 

We will now use this result to prove the second theorem in this section. 

Proof of Theorem 3.2. It follows from Lemma 3.7 and Theorem 3.1 that all but finitely 
many limit points of the zeros of the sections sn[F](Nz) with $t(z) ^ must lie on the 
curve D a ,b as defined in part (i) of Theorem 3.2. That every point on D a j, is such a limit 
point follows from the choice of the subsequence {-/V} which ensures that the sequence of 
sections {sjv[-F](,z)} has a positive fraction of zeros in any sector with vertex at the origin. 
It is straightforward to show that for any < 9 < 2tt there is a unique r > such that 
re 10 € D a ,b- This proves part (i) of Theorem 3.2. 

From the asymptotic expansion for F(nz) in Lemma 3.4 we see that the limit points of 
the zeros of F(nz) must lie on the imaginary axis. Further, we claim that F has infinitely 
many zeros, so every point on the imaginary axis is such a limit point. Indeed, if F had 
finitely many zeros then, since F has order 1, the Hadamard factorization theorem would 
tell us that we can find a constant a and a polynomial p such that 

(9) F(z) = e az p(z). 

But if y e K then 

\F(W)\ -- 



b 

< I \<p(t)\dt. 



b 

ip(t)e iyt dt 

-a 

Taking the limit as z — > ioo in (9) we see that we must have 9(a) > 0, and taking the limit 
asz-> —ioo we must have 3(a) < 0, a contradiction. Thus the claim is verified. 

If {zn} is a sequence of complex numbers such that sn[F](Nzn) = for all N which has 
a limit point on {z S C : di(z) = and 9(z) > l/(ec)} then by equation (6) and Lemma 3.6 
we must have F(Nzn) oo. We just showed that F(iny) is bounded, so such a sequence of 
zeros cannot exist. Hence any limit points on the imaginary axis must satisfy Q(z) < l/(ec). 

Further, if {zjv} is a sequence of zeros which has a limit point in \z\ < l/(ec), then by 
equation (6) and Lemma 3.6 we must have F(Nzn) — > 0. In other words, the zeros {zn} 
must approximate the zeros of F(Nz). Thus part (ii) is proved. 

The truth of part (iii) follows from combining the facts that a) all but finitely many limit 
points with 5R(z) ^ must lie on D a b and b) any limit points with 5R(z) = must lie on 
Amag- This completes the proof of Theorem 3.2. □ 



4. Special cases of the exponential integrals 

We mentioned in the introduction that the confluent hypergeometric functions with b > 1 
studied by Norfolk [11] are a special case of the exponential integrals studied in this paper. 
Indeed, when b > 1 we have the integral represenation 

iFi(l;6;z) = {b - 1) f (1 - t) b ~ 2 e zt dt. 
Jo 

Our result extends some of Norfolk's results to the case of complex 6 with $ft(6) > 1. 
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Another important special case of the exponential integrals is the class of Bessel functions 
of the first kind, defined by 



:\"x— (-l) fc fz^ 2k 



2/ ^T{k + l)T{k + a + l) \2 
For $t{a) > —1/2 we have Poisson's integral representation (see, e.g., [17] 

Uz) = -, ffi" m t (l-t 2 r 1/2 e tzt dt, 



r(«+|)r(|) 

which is of the required form once z is replaced with — iz. As far as we can tell, asymptotics 
for the zeros of sections of the Bessel functions have not been previously studied. We 
state the result formally as a corollary to Theorem 3.1 and Theorem 3.2. This corollary is 
illustrated in Figure 3. 

Corollary 4.1. Let J a be the Bessel function of order a with $t(ct) > —1/2, and for n even 
let 

^/2 



AJ a ]{z) - ( 2 ) E r(fc + 1 ) ( r (fc + a + i) ( 2 ) 



be its n th section. Then, for a subsequence {N} of the indices {n} as defined in Theorem 2.1, 
all but finitely many limit points of the zeros of the normalized sections s^r[J a ](Nz) lie on 

D(J) = {z e C : 3(z) > 0, |*| < 1, and \ze 1+lz \ = 1} 

U {z e C : < 0, \z\ < 1, and \ze^ lz \ = 1} 

U{i£l: -l/e<x< 1/e} . 

Conversely, every point of D(J) is a limit point of zeros. Zeros of sn[Jo\(Nz) in > 
which do not converge to the point z = i satisfy 

\^\=l + l ^ + 0(l/N) 

as N — > 00, and zeros in ^s(z) < which do not converge to z = —i satisfy 

1-1 IoeTV 

l ze l = 1 + ^\r + 0(1/A ° 

as N — > 00 . Furthermore, if a is real, the modified indices {N} in this result may be replaced 
everywhere by the original even indices {n} and, in this case, all limit points of the zeros of 
the normalized sections s n [J a ](nz) lie on D(J). 

Proof. The only part of Corollary 4.1 which does not follow immediately from 
Theorem 3.1 or Theorem 3.2 is the claim that if a is real, the modified indices {N} may be 
replaced by the original even indices {n} and, in this case, all limit points of the zeros of 
the normalized sections s n [J a ]{nz) lie on D{J). 
To prove this claim we write 



s n [J a ](-inz) = P n (z 2 ) 



where 

n/2 



o 4 fc r(fc + i)r(fc + a + 1) 



z k . 
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Figure 3. LEFT: Zeros of the normalized section s 80 [Jo](80z) and the 
limit curve D{J) as defined in Corollary 4.1 . 

By the Enestrom-Kakeya theorem (see, e.g., [9, ch. 7] and [1]), all zeros of P n {z) satisfy 

, 2n + 4 r(^+a + 2) 

\z\ < • — — — 

1 1 " n 2 r(f + a + l)' 

We apply Stirling's formula for the Gamma function to see that 

r(f +a + 2) n 
r(f +a+l) ~ 2' 

whence 

2n + 4 r(f +a + 2) 
n 2 ' r(f +a + l) ~^ 

as n — ¥ oo. This tells us that the limit points of the zeros of the polynomials P n , and hence 
of the sections s n [J a ](nz), lie in the closed unit disk. This replaces the use of Theorem 2.1 
and hence Lemma 3.7 in the proofs of Theorem 3.1 and Theorem 3.2. Indeed, Theorem 2.1 
is the origin of the restriction of the indices to the subsequences {^V}. This completes the 
proof of the corollary. □ 
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5. Discussion and future work 

The exponential integrals in this paper were studied in part because of their simple 
definitions. It should be noted that their integral form facilitated symbolic manipulations 
and allowed for the use of standard asymptotic techniques. But they were also studied 
because their definition was flexible enough to allow the zeros of their sections to display 
some unique behavior, perhaps the most notable of which is the tendency for the zeros to 
approach the limit curve from the interior or the exterior depending on the orders of the 
singularities/zeros of the integrand tp at the endpoints of integration (see Remark 3.3). 

As mentioned in Section 1, this paper is the first in a series of two. In the forthcoming 
second paper we will complement Theorem 3.1 by studying the asymptotics of the zeros 
which approach the points of the limit curve D a jy on the circle \z\ = 1/c, which occur at 
z = 1 / b or z = —1 / a (or both). In general, the zeros approach these points more slowly than 
they approach the body of the curve. In fact, the unique behavior of the zeros described in 
the previous paragraph will play an important role in the analysis. Just as in this paper, 
results about the zeros of sections of the Bessel functions will appear as special cases. 
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